A transient solution for droplet deformation under electric fields 
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A transient analysis to quantify droplet deformation under DC electric fields is presented. The full 
Taylor-Melcher leaky dielectric model is employed where the charge relaxation time is considered 
to be finite. The droplet is assumed to be spheroidal in shape for all times. The main result is 
an ODE governing the evolution of the droplet aspect ratio. The model is validated by extensively 
comparing predicted deformation with both previous theoretical and numerical studies, and with ex- 
perimental data. Furthermore, the effects of parameters and stresses on deformation characteristics 
are systematically analyzed taking advantage of the explicit formulae on their contributions. The 
theoretical framework can be extended to study similar problems, e.g., vesicle electrodeformation 
and relaxation. 



I. INTRODUCTION 

When a liquid droplet suspended in another immis- 
cible fluid is subject to an applied electric field, it un- 
dergoes deformation due to the electrostatic stresses ex- 
erted on the interface. Extensive research on this phe- 
nomenon has been conducted to study the deformation 
due to its relevance in a variety of industrial applications, 
including electrohydrodynamic atomization electro- 
hydrodynamic emulsification [H, and ink-jet printing [sf, 
among others. Historically, the deformation dynamics is 
divided into two regimes: electrohydrostatics (EHS) and 
electrohydrodynamics (EHD). In the first, EHS deforma- 
tion, the droplet is idealized as a perfect conductor im- 
mersed in a perfect insulating fluid; or both of the fluids 
are treated as perfect dielectrics with no free charge 0- 
1^. For this case, the electric field only induces a normal 
electrostatic stress, which is balanced by surface tension, 
and the final equilibrium shape is always prolate. At the 
steady state, the hydrodynamic flow is usually absent. 
In the second, EHD deformation, both fluids are consid- 
ered to be leaky dielectrics [tJ, Il0l420| . For this case, 
when an electric field is applied, free charges accumulate 
on the droplet surface which induces a tangential elec- 
trostatic stress in addition to the normal one. Driven 
by this force, the fluids inside and outside the droplet 
present toroidal circulations and a viscous stress is gen- 
erated in response to balance the tangential electrostatic 
stress [l^. The droplet deforms into either a prolate or 
an oblate spheroid shape depending on the specific elec- 
trical properties of the fluids. With different electrical 
properties, the effects of the electrostatic and hydrody- 
namic stresses on droplet deformation are distinctive. 

This work focuses on a solution method for problems of 
the second kind, namely, EHD deformation. This type of 
problem is more challenging to solve. In the literature, all 
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theoretical solutions were obtained largely under two spe- 
cific assumptions: (i) The deformations are small. The 
analysis is performed by assuming that the equilibrium 
shape of the droplet only slightly deviates from spheric- 
ity. Solutions using this assumption can be found in 

[1- 

Il2| . (ii) For large deformations, the shape is assumed 
to be spheroidal during the entire deformation process. 
Results using this assumption are given in Ts]. When 
compared with experimental data, predictions from the 
small-deformation theories always quantitatively under- 
predict the aspect ratio especially when the deformation 
is large. In contrast, the large-deformation theory has a 
better agreement both qualitatively and quantitatively. 
In all of the above, the theoretical analysis only leads 
to solutions in the steady state. The Taylor-Melcher 
leaky dielectric model |2l| - [23| with the assumption of 
instantaneous charge relaxation has always been used. 
On the other hand, the theoretical analysis of transient 
droplet deformation seems to attract less attention. Only 
Dubash and Mestel [9| developed a transient deformation 
theory for an inviscid, conducting droplet. This analysis, 
which solves a EHS deformation problem, is not applica- 
ble to study EHD deformations. In general, to fully solve 
the transient EHD problem, numerical simulations have 
been employed [H-lilliO]. 

In this work, we present a transient analysis of droplet 
deformation under direct-current (DC) electric fields. 
Following Bentenitis and Krause [l^, we assume the 
droplet remains spheroidal in shape. The full Taylor- 
Melcher leaky dielectric model is employed where the 
charge relaxation time is considered finite. In this frame- 
work, instantaneous charge relaxation is treated as a spe- 
cial limiting case. This generalization allows direct com- 
parison with experimental data which were usually ob- 
tained in fluids with very low conductivities [8]. The 
main result is an ordinary differential equation (ODE) 
governing the evolution of the droplet aspect ratio. The 
availability of this equation allows us to explicitly analyze 
the effects of parameters and stresses on the deformation 
characteristics. The model is validated by extensively 
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comparing predicted deformation with both previous the- 
oretical and numerical studies, and with experimental 
data. 



II. THEORY 

A schematic of the problem configuration is shown 
in Fig. [IJa). An uncharged, neutrally buoyant liquid 
droplet of radius tq is suspended in another fluid, and is 
subject to an applied electric field of strength Eg. We 
assume that the fluids are immiscible leaky dielectrics 
with constant electrical and mechanical properties, tr, e, 
and fj, are the electrical conductivity, permittivity, and 
fluid viscosity, and the subscripts i and e denote inter- 
nal and external, respectively. Under the influence of an 
applied electric fleld, free charges accumulate at the inter- 
face, which induces droplet deformation and EHD flows 
both inside and outside the droplet. Taylor ^lOj predicted 
that droplets may deform into prolate or oblate shapes 
depending on the electrical properties of the fluids. In 
the following analysis, we focus on developing a solution 
for prolate deformations, whereas a solution for oblate 
deformations can be pursued in a similar manner (not 
presented here). 

We assume that the droplet remains spheroidal in 
shape throughout the process. This approximation is 
consistent with experimental observations by Ha and 
Yang and Bentenitis and Krause . Following Tayor 
a|, Bentenitis and Krause [18], and Dubash and Mestel 
9|, the natural coordinate system to analyze this prob- 
lem is the prolate spheroidal coordinate system, and a 
schematic is shown in Fig. [Ub). The geometry is as- 
sumed to be axisymmetric about the z axis, which aligns 
with the direction of the applied electric field. The 
spheroidal coordinates -q) are related to the cylindrical 
coordinates (r, z) through the equations: 



(a) 

Eq 




ri = const 



(1) 



FIG. 1. (a) A schematic of the problem configuration, (b) 
The prolate spheroidal coordinate system. 



with deformation. The critical idea of the current anal- 
ysis is to express all variables, e.g., the electric potential 
and the stream function in terms of ^o- 

In what follows, we will solve the electrical problem 
first, followed by a solution of the hydrodynamic problem. 
An ODE for is obtained by applying both the stress 
matching and kinematic conditions. 



r = c^ie-m-v')- (2) 

Here c = Va? — b'^ is chosen to be the semi-focal length 
of the spheroidal droplet, and a and b are the major and 
minor semi-axis, respectively. The contours for constant 
^ are spheroids, and ^ S [1, -f-oo). The contours for 
constant r] are hyperboloids, and rj £ [—1,1]. The surface 
of the prolate spheroid is conveniently given as 

e = Co^-. (3) 

c 

For the derivation below, we further assume that the vol- 
ume of the droplet is conserved. We subsequently obtain 

a = ^o(l-^o"')"^ 6 = ro(l-Ccr')^- (4) 

Therefore, the droplet geometry is completely character- 
ized by a single parameter, which evolves in time along 



A. The electrical problem 

The electric potentials inside and outside the droplet 
obey the Laplace equation according to the Ohmic law 
of current conservation with uniform electrical conduc- 
tivity: 

V^cj), = V^cj)e = 0. (5) 
The matching conditions at the interface are 

||V./)-t|| =0, at^^fo, (6) 

||-||aV0-n||=O, ate^Co. (7) 

Here 1 1 • 1 1 denotes a jump across an interface, and t and 
n are the unit tangential and normal interfacial vector. 
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respectively. q= ||—eV0-n| | is the surface charge density. 
Note that in Eq. ([7]), we have included the displacement 
current, dq/dt. This term is particularly important for 
fluids with very low conductivities (for example, those 
used in Ref. 8]) such that the interfacial charging time 
becomes comparable to the deformation time. However, 
we have ignored the effect of surface charge convection 
which is shown to be small by numerical simulations [l6| . 
Equation ([7]) can be rewritten in terms of the electric 
potentials. 



Ope _ £i o<Pt 
c dS, c J dt 



_e^a0e ei_9^\ dc 

h^c dS, h^c J dt 

ate = ^o- (8) 



Here is a metric coefficient of the prolate spheroidal 
coordinate system. The displacement current consists of 
two parts, represented by the first three terms on the 
LHS of the above equation. The first two terms result 
from a change in the droplet shape, and the third re- 
sults from the charging process as if the shape remains 
unchanged. Here, we assume that the first term is negli- 
gible compared with the other two, and Eq. (|8]) can be 
further simplified to be 



92 



dc 
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dS^dt c d( dt 'd^dt 









-0, 








Tl 


Qi(eoKo i' 


dp 




T2 









Ci d4ii dc 

at e - Co- (9) 



Indeed a consistency check a posteriori justifies the sim- 
plification. Far away from the droplet surface the electric 
field is uniform, 



V0e Eqz, at ^ 



oo. 



(10) 



We also require that (j>i remains finite at ^ = 1. For the 
initial condition, we assume both the electric potential 
and the normal component of the displacement vector 
are continuous: 



or 



ate = eo, t = 0. (11) 



Solutions for the electric potentials have been obtained 
previously without including the displacement current 
|10l [l8| . With its inclusion the approach is similar and 
the results are 



(l)e = Eoro [-Ae + aQi(?)]r/, 



= Eorol3^V- 



(12) 



(13) 



Here, Qi{^) is a Ist-degree Legendre polynomial of the 
second kind. A = c/tq is the dimensionless semi-focal 
length. The coefficients a and (3 are determined by the 
interfacial matching conditions ^ and © which gives 



a = 



Qi(eo) 



:(eo)gi(eo)-Q?(eo) 



Qi(eo)eo 



1 

(Jr 



/3 + A 



Q'i(Co)go 



T2 



Qi(^o) 




P^o + ACo 

Qi(So) 



Oi(eo)eo 



(14) 



+ 



Qi(eo)A 
i(go)Qi(eo)^Q?(^o) 

Q?(?o) 



dX d<^o 
d^o dr 

A^ = 0,(15) 



a(0) = A^o (er - 1) , /3(0) 



A (Qi(Co) - (9'i(Co)fo) Eq. (dll) can be simplified to be 



erQi(eo)eo-Qi(eo) 



(16) 

Here = ee/e^ and Ur = <Je/<Ji are the permittivity ratio 
and the conductivity ratio, respectively, ri = Ce/cTe is an 
electrical charging time. T2 = r^^e/l ^ characteristic 
flow timescale used below in the hydrodynamic problem, 
and 7 is the coefficient of surface tension. In the above 
equations, a dimensionless time t = t/T2 has been used. 
In general, Eq. (|15p needs to be integrated together with 
an ODE for to obtain a and /3. However, in the limit 
of instantaneous-charge-relaxation time, ti/t2 — > 0, and 



Qi(Co)^o 

Qi(eo) 



1 

(Jr 



/3 + A 



Qi(eo)^o 

Qi(eo) 



- 1 



= 0. 



(17) 



This result is equivalent to a solution employing the sim- 
plified boundary condition | |(tV0 • n| | = in place of Eq. 
©• 

The normal and tangential electrostatic stresses are 
given by. 



Ei 



S^j-i — eE^E,j, 



(18) 



where E^ — —{d(t)/d^)/h^ and Ejj = —{d(j)/dri)/hrj are 
the normal and tangential electric fields, respectively, 
is a metric coefficient of the prolate spheroidal coordi- 
nate system. These stresses can be evaluated with the 
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solutions p2|) and p3|) . and will be used in the stress 
matching conditions below. 



B. The hydrodynamic problem 

In the regime of low-Reynolds-number flow, the gov- 
erning equation for the hydrodynamic problem can be 
rewritten in terms of the stream function, ip, as 



E^iP = 0. 



(19) 



Here, the expression for the operator can be found in 
Dubash and Mestel 19] and Bentenitis and Krause [l8| . 
The stream function is related to the velocity components 
as 



1 dtfj 



1 dtp 
hr/he dr] 



(20) 



he is a metric coefficient of the prolate spheroidal coor- 
dinate system. At the interface, u and v represent the 
tangential and normal velocities, respectively, and they 
are required to be continuous 



at ^ = Co- 



(21) 



In addition, we prescribe a kinematic condition relating 
the interfacial displacement to the normal velocity. 



vis. = Co, v) 



(1 - 



-5/6 



1 - 3r;2) 



dt 



(22) 



The total force on the interface resulting from the elec- 
trical stress, the hydrodynamic stress, and the surface 
tension should be balanced at every point. However, 
this constraint is impossible to satisfy exactly within 
the framework of spheroidal deformation. Various au- 
thors develo ped reduced stress-balance conditions in- 
stead 0, [glllSl . Here we follow the integrated formulae 
proposed by Sherwood and Dubash and Mestel |[9|] 



in 



(23) 



ds = 0. 
(24) 

Equations ([23l) and (f24| represent a global balance of 
the tangential and normal stresses, respectively derived 
from energy principles. Here T denotes the hydrody- 
namic stress, Ri and i?2 are the two principal radii of 
the curvature, and the integration is carried over the in- 
terface. 

The general solution to was proposed by Dassios 
et al. |24| using the method of semi-separation: 

OG 

+ J2^9n{OGn{v) + KiOH^m ■ (25) 

n=2 



Here G„ and i?„ are Gegenbauer functions of the first 
and second kind, respectively. f/„ and /i„ are linear com- 
binations of Gn and The detailed exp ressions for 
Gn, Hn, gn, and /i„ are found in Rcf. [24!|. Interested 
readers are referred to Ref. for further details. After 
considering that the far field is quiescent, and that the 
velocities remain finite at C = 1, the stream functions can 
be simplified to be 

oo 



n=l 



oo 



+Al:^XlH2n+3iO] G2n+liv), (26) 



^» = 1^ [Bi:^iG2n-iio + Bi:iiG2n+iio 

n=l 

+ Bl:+lG2n+3{0] G2n+liv), (27) 

where A and B are unknown coefficients satisfying the 



relations Al^+l = AllX\, BllH 



i?2tT3 . In general, 
these coefficients are inter-dependent, and the full solu- 
tion can be obtained only with the entire infinite series. 
Here we seek a truncated solution as an approximation, 

tP,= [A\H,{^) + AlH:i{0]GM, (28) 



= [i?|G3(0 + BlG^iO] GM- (29) 

Indeed, G^^r]) gives a functional form in rj confirming 
with that in Eq. (P^ . which can be rewritten as 



2c'V^^ro (1 - SoT'^' G'M d^o 



v{S ^ Co, v) 



3Co^ 



hfjhg dt 



(30) 

This agreement in part validates the spheroidal shape 
assumption: the shape represents the leading mode in 
the infinite series. 

Equations (|21j|24p are combined to solve for the five un- 
known variables, namely, A^, A^, JSf, Bf, and Co- Specif- 
ically, Eqs. (|2T|) and (|22l) are first used to eliminate the 
Ah Bl Bl 



Al - H3{^o)Al ~ M 



rfCo 
dt ' 



„3 -g5(Co)g;(Co)Ai+G;(Co)Mf^ 

^3- ^ 



„5 G3(Co)g;(Co)^g-G:3(Co)Mf^ 

N ' 



(31) 



(32) 



(33) 



where M = 2r3/3(Co' - Co), and N = G3(Co)G;(Co) 
G3(Co)G5(Co)- Further considering Eq. 
press in terms of Co, 



we can ex- 



{^0/3' - er(A - ag;(^o))(ACo - aQi(Co))} hi{(,o) - {(m. - l)/i2(eo) + /i3(Co)} 



-Mi {Air /l4 (Co) + /l5(Co)} 



(34) 



I 

where //^ = fJ-e/fJ-i is the viscosity ratio. The detailed This expression is inserted into Eq. (j24p to obtain the 
expressions of /ii(Co)~/i5('?o) ai'G found in the Appendix. final result, an ODE governing the evolution of the ^O; 



dr F 



D f fC \ , D Mr/22(Co) + /23(Co) r /c \ 
Mr/14 (4o) + /l5(to) 



(35a) 



Q 



N 



(A - aQi(eo))' + (A - agi(Co)/Co)' - 2/3Ver 



(35b) 



(A - agi(Co))(A - agi(eo)/So) - /^V^r 



(35c) 



The detailed expressions of /2i(Co) ^ /24('?o), and 
are also found in the Appendix. The coefficients a 
and /? are given by Eqs. ^(T^ and p^ . respectively. 
Cqe = is the electric capillary number. In 

Eq. p5ap , the three terms in the numerator on the RHS 
represent the contributions from the normal stress, the 
tangential stress, and the surface tension, respectively. 
At the equilibrium, the balance of the three forces deter- 
mines the final shape. The leading coefficients Qn and 
Qt arise exclusively from the electrostatic stresses, and 
can be used to estimate their respective influence on de- 
formation. In the limit of instantaneous relaxation, and 
by considering Eqs. (fT4|) and ([TTl) . Qn and Qt can be 
simplified to be 



CaEK'^{'7l+l-2al/er).QT = CaEK^(Tr{l-ar/er), 

(36) 



K I 



?l(Co) - CTr^oQ'li^o)' 



(37) 



For this case, the evolution of is governed by a single 
timescale, T2. Once is obtained by solving the Eqs. 
([T5|) and p5ap . the aspect ratio is calculated by the for- 
mula 



i^(l-eo- 



(38) 



III. COMPARISON WITH PREVIOUS RESULTS 

In this section, we compare our model prediction exten- 
sively with results from previous work. The comparisons 
with theoretical/numerical results and experimental data 
are respectively presented in Sees. IIII Al and IIIIBl 



Comparison with previous theories and 
simulation 



We first consider the equilibrium shape, and com par e 
our results with those from Bentenitis and Krause |l8| . 
For this case, the LHS of Eq. p5a|) is simply set to zero, 
resulting in the so called discriminating equation. 



QNf2l{^0) + QT 



Mr/22(Co) + f 23(^0) 
Mr/l4(Co) + fl5{£,()) 



/24(eo). (39) 



Here Q^ and Qt are given by Eq. ([551) . ^0 is solved 
as a root(s) of this equation from which the equilibrium 
aspect ratio, a/b, can be obtained. Equation (j39p shows 
that the equilibrium shape is only determined by the di- 
mensionless parameters Cag, ar, e^, and ^r- A com- 
parison with the theoretical prediction by Bentenitis and 
Krause [3 is shown in Fig. [21 Note that in this earlier 
work, the authors solved for the equilibrium shape di- 
rectly without obtaining the transient solution. A good 
agreement is observed, although a different stress match- 
ing condition has been used by Bentenitis and Krause 
13 [see Eqs. (38) and (45) therein]. 

We next compare with the results from Dubash and 
Mestel Q. In this work, the authors developed a theo- 
retical model, also with the spheroidal shape assumption, 
to predict the transient deformation of a conducting, in- 
viscid droplet immersed in a viscous, nonconductive solu- 
tion. This special consideration leads to significant sim- 
plifications: both the electric and hydrodynamic fields 
are absent within the droplet. In addition, at the equilib- 
rium state (if one is permitted), the hydrodynamic flow 
outside the droplet is also quiescent, giving rise to the 
phenomenon termed EHS. 

In our generalized framework, the solution for this case 
is simply achieved by setting ar ^ and /i^ — > c» in 
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FIG. 2. The equilibrium aspect ratio as a function of electric 
capillary number. The parameters are Or ~ 1.19 x 10"'^, 
e,. = 3.24 X 10-\ and = 7.33 x 10"^ 

Eqs. (I35ap and Note that ^ directly leads 

to instantaneous charge relaxation. The resulting com- 
parisons are shown in Fig. |3] in which the aspect ratio 
(a/&) is plotted as a function of time for four different 
electric capillary numbers {CaE)- For the two lower val- 
ues of CaEi the current model has excellent agreement 
with both the theoretical and numerical predictions by 
Dubash and Mestel [9| [Fig. [Sta)]. For these Gas val- 
ues, final equilibria are achieved. As Cas increases [Fig. 
[3]Jb)], the deformation becomes unstable and an equilib- 
rium shape is no longer possible. The rapid expansion 
with a sharp slope at the later stage preludes droplet 
breakup. For these two cases, the theoretical models still 
agree with each other, whereas some discrepancies exist 
with respect to the numerical simulation, in particular 
for Ca£;=0.206. However, this discrepancy is in general 
only noticeable when the CaE number is above and very 
close to the critical threshold of breakup (CaE ~0.2044 
for the case studied), due to a slight underprediction of 
the rate of deformation by the theoretical models. A sim- 
ilar trend is observed when comparing with the numerical 
simulation by Hirata e< al. jl7j (not shown). Overall, our 
model can serve as a good approximation to the numer- 
ical model which is considered more accurate. 



B. Comparison with experimental data 

The main source of experimental data comes from Ha 
and Yang Q. We also begin with an examination of 
the final aspect ratio when an equilibrium shape can 
be achieved. Figure H] shows the equilibrium aspect ra- 
tio of a castor oil droplet immersed in silicone oil from 
Ha and Yang as well as predicted by various mod- 
els. The current prediction is shown as a solid line, 
whereas the results from first-order (To| and second- 
order [l2l| theories are shown as dot-dashed and dashed 



(a) 1.8 




20 40 60 80 100 

r 




20 40 60 80 100 



T 

FIG. 3. The deformation of a conducting droplet in a highly 
viscous medium, (a) Cas = 0.18 and 0.204. (b) Cas = 0.206 
and 0.21. The dimensionless time r is defined as r = i/T2, 
where T2 = ro/ie/7. 



lines, respectively. Following Lac and Homsy [20|, we 
rescale CaE to best match Ajayi's second-order correc- 
tion. This rescaling is equivalent to adjusting the sur- 
face tension from 7 = 3.3 x 10~^ N/m used by Ha and 
Yang 8] (which is a fitting parameter in that work) to 
7 — 4.3 X 10^'' N/m. The latter value is close to the lower 
bound, 7 = 4.5 x 10"'^ N/m, measured by Salipante and 
Vlahovska In addition, we use Ur = 0.03 accord- 

ing to the measurements by Torza et al. |11|. Vizika and 
Saville [l^, and Salipante and Vlahovska |25j|. which is 
slightly different from the value of = 0.04 used by 
Lac and Homsy [2^. The results show good agreement 
between the current model and the experimental data. 
Most importantly, our theory correctly predicts a criti- 
cal CaE (~0.244) for droplet breakup. In contrast, the 
small deformation theories can not capture this critical 
phenomenon. 

We have also compared our theoretical prediction with 
the experimental data from Bentenitis and Krause fisj . 
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FIG. 4. The prediction from current model is compared with 
the small deformation theories [lol.[l^ and experimental data 
[^. The parameters are ar = 0.03, tr = 0.73, and /x^ = 1.14. 



which measured the equilibrium aspect ratio of a DGEBA 
droplet immersed in a PDMS solution. Since our result is 
in good agreement with the theoretical prediction in the 
same work (see Fig. [5]), which in turn agrees well with 
the data, the comparison is not shown here for brevity. 

Next, we will compare the transient solution from our 
model with data from Ha and Yang Q. In Fig. [SJa), 
the data is extracted from Fig. 3 in the latter work, 
which captures the deformation of a water droplet in sil- 
icone oil. The droplet is fitted with an ellipse at every 
instant, based on which the aspect ratio is calculated. 
A 10% fitting error is estimated, and is shown as er- 
ror bars in Fig. El^a) [the same approach is adopted to 
extract the data presented in Figs. El^b) and[B]. The 
model prediction is calculated with Eqs. pSal) and (l36t . 
and with cr^ = 1 x 10^^ = 3.55 x 10~^ = 1000, 
Eq = 3.2 kV/cm, ro = 0.25 cm, and = 0.98 Pa • s 
all directly taken from Ha and Yang [8]. For medium 
permittivity, we use te — 2.478 x 10~^^ F/m following 
the measurements by Torza et at. [lll |. Vizika and Sav- 
ille [l3|, and Salipante and Vlahovska [1^. For surface 
tension, we use 7 = 3.037 x 10~^ N/m, which is consis- 
tent with the values reported by Torza et al. [ll[ and 
Vizika and Saville [l^. In this case, the model is able to 
predict the deformation process with good quantitative 
accuracy. In Fig. [SJb), a similar comparison is shown for 
a water-ethanol droplet in silicone oil. The data is based 
on Fig. 4 in Ref. Q. For our calculation, ar — I x 10~^, 
Cr = 0.05, = 23.3, Eo = 4.5 kV/cm, ro = 0.14 cm, 
fie = 0.98 Pa • s, and = 2.478 x 10"" F/m. Because 
the droplet is doped with polyvinylpyrrolidone (a poly- 
mer solution), the surface tension is not directly avail- 
able, and is used as a fitting parameter instead to gen- 
erate the best agreement between theory and data. The 
resulting value is 7 = 3.432 x 10~^ N/m, 11% higher than 
that for water/silicone oil which is used in Fig. [SJa). 

In contrast to the regime of instantaneous charge re- 
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FIG. 5. Comparison of transient droplet deformation, (a) 
A water droplet in silicone oil. The parameters are a,- = 

= 1000, Eo = 3.2 kV/cm, 
ee = 2.478 x 10~" F/m, 



1 X 10"^ er = 3.55 X 10-^ ^Ir 
ro = 0.25 cm, fie = 0.98 Pa ■ s 



and 7 = 3.037 x 10" N/m. (b) A water-ethanol droplet in 
silicone oil. The parameters are (Tr = 1 x 10~^, er = 0.05, 
fir = 23.3, Eo = 4.5 kV/cm, ro = 0.14 cm, fie = 0.98 Pa • s, 
= 2.478 X 10"" F/m, and 7 = 3.432 x 10"^ N/m. 



laxation examined in Fig. [5l Fig. [B] represents droplet 
deformation in the finite-charging-time regime. The data 
is extracted from Fig. 7 in Ref. [8]. For this case, the 
droplet is made of castor oil, and is immersed in sili- 
cone oil. The extremely low conductivities of these me- 
dia lead to a charging time (^seconds) comparable to 
the deformation time, and the full model, Eqs. (|35ap - 
pScp . has to be used. For our calculation, ar — 0.03, 
er = 0.73, fir = 1-14, Eo = 3.2 kV/cm, ro = 0.16 cm, 
fie = 0.9 Pa-s, Ee = 2.478 x 10"" F/m, and 7 = 
5 X 10""^ N/m. Note that the values for the surface 
tension and the conductivity ratio follow the measure- 

Vizika and Saville and 



ments by Torza et al. [11 



Salipante and Vlahovska [25| which are believed to be 
more accurate than the original values of ar =0.1 and 
7 = 3.3 X 10~^ N/m given by Ha and Yang In ad- 
dition, the actual conductivity of silicone oil varies from 
10-1° S/m to 10"" S/m in the literature (H-H^. In Fig. 
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FIG. 6. Droplet deformation in the limit of extremely low conductivities. The parameters are cr^ — 0.03, tr = 0.73, Hr = 1.14, 
Eo = 3.2 kV/cm, ro = 0.16 cm, fie = 0.9 Pa • s, te = 2.478 x 10"" F/m, and 7 = 5 x 10"^ N/m. The best agreement between 
the data and the theory is found for CTo = 7 x 10~^^ S/m. For reference, the dotted line shows the calculation according to the 
instantaneous-charge-relaxation (ICR) model. 



[SI we show the calculation with three representative val- 
ues within this range, namely, CTg = 1 x 10~^°, 7 x 10~^^, 
and 1 X 10""'^^ S/m. The best agreement is found for 
ae — 7 X 10^ ^■^ S/m. For comparison, the calculation 
according to the instantaneous-charge-relaxation model 
[Eqs. (|35ap and ([36]) ] is also shown, and is denoted by 
ICR. This simplified model clearly overpredicts deforma- 
tion by a significant degree. 

In general, our model agrees well with experimental 
data in both steady and transient states, and for a large 
parametric range. These comparisons provide a strong 
validation for our model. 



IV. THE EFFECTS OF STRESSES ON 
DEFORMATION 

In this section, we demonstrate the utility of our theo- 
retical results by analyzing in-depth the governing equa- 
tion. For simplicity, we focus on the regime of instan- 
taneous relaxation, where and Qt are given by Eq. 
P6p . A main contribution of the current work is that 
Eq. (|35aP clearly separates the effects by different forces. 
In the numerator of the RHS, the three terms represent 
respectively the effects of the normal stresses (both elec- 
trical and hydrodynamic) , the tangential stresses (both 
electrical and hydrodynamic), and the surface tension. 
Furthermore, all the functions in this equation are posi- 
tive (/i4, /i5, /21 — /24, F), such that the signs of Qat and 
Qt completely determine whether the normal and tan- 
gential stresses would promote or suppress deformation. 
Due to the inverse relationship between ^0 and the aspect 
ratio, a/b [see Eq. (|38p ]. a positive Qn or Qt indicates a 
positive contribution. Evidently, surface tension always 



resists deformation. Because Qn and Qt depend exclu- 
sively on the electrical properties in a simple manner [see 
Eq. p6p ]. their influences can be conveniently analyzed 
using a phase diagram shown in Fig. [T] The dashed 
and dotted lines correspond to Q^v = and Qt = 0, 
respectively. These lines separate the phase space into 
three regimes, where N and T denote the normal and 
tangential stresses, and the superscripts '+' and '— ' de- 
note a positive or negative contribution to deformation, 
respectively. In addition, the solid line is obtained by 
solving for the root of Taylor's discriminating function 
p^ . which separates the prolate (denoted by 'Pr') and 
oblate (denoted by 'Ob') regimes [this line can be equiv- 
alently obtained by looking for the steady-state solution 
of a/6 = 1 from Eq. ([55)) ]. 

Figure [7| can be used to shed light on the physical pro- 
cesses governing deformation. First, the line for Qt = 
separates the T"*" and T~ regimes, which corroborates 
with the previous results [20|. On this dividing line, 
the velocity field becomes zero, so does the tangential 
electrical stress. In Ref. [13], the viscosity ratio has op- 
posite effects on deformation in the r+ and T" regimes. 
This behavior is clearly explained by Eq. (|35ap . Second, 
there is a small region within the oblate regime, namely, 
the area between the solid and dashed lines where Qat is 
positive. This suggests that the normal stress still tends 
to stretch the droplet along the direction of the applied 
field. However, because Qt is negative, the tangential 
stresses overcome the normal stresses, and stretch the 
droplet into an oblate shape. This new insight is not 
available from previous analysis or simulations. 

Third, in the prolate regime where Qn is always pos- 
itive, the sign of Qt leads to different deformation be- 
havior. Figure [5| shows the equilibrium aspect ratio as a 
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FIG. 7. Phase diagram for droplet deformation. Here the 
dotted hne is calculated by satisfying Qt = 0. The solid line 
is calculated by solving for the root of Taylor's discriminating 
function. The dashed line represents Qiv = 0. Pr and Ob 
denote prolate and oblate deformation, respectively. A'^ and 
T denote the effect of normal and tangential stresses, respec- 
tively, and a '+' or '— ' sign denotes facilitating or suppressing, 
respectively. 

function of Cas for three specific cases. Note that the 
new variable 



In this new definition, D — corresponds to a/b = 1, 
and D = 1 corresponds to a/b ^ oo. For all three cases, 
er = 10 and fj.r = 1- For ar = 0.05, Qt > 0. We ob- 
serve hysteresis, and D approaches 1 rapidly in the upper 
brunch. The cases of = 1 and ar = 30 correspond to 
Qt = and Qt < 0, respectively. In general, as Qt de- 
creases, the deformation becomes weaker for comparable 
CaE values. Most interestingly, for ar — 30 {Qt < 0), D 
converges to a value less than 1 in the limit of Ca^ -> oo. 
This means that even for the very large applied electric 
field strength, a finite equilibrium aspect ratio can be 
achieved. We emphasize this scenario is only possible 
in the T~ regime. For large Eq values, corresponding 
to large Cue, the resistive effect from surface tension is 
negligible, and the only way to obtain a finite equilib- 
rium aspect ratio is therefore by balancing the normal 
and tangential stresses. Since Qn positive, Qt has to 
be negative. 




0.2 0.4 0.6 0.8 1 



FIG. 8. The behavior of equilibrium droplet deformation in 
different regimes. For ar = 0.05, Qt > 0; Ur = 1, Qt = 0; 
(7r = 30, Qt < 0. As Cue — ^ oo, an equilibrium shape is only 
possible in the T~ regime. Other parameters are tr = 10 and 
Mr = 1. 



is assumed to be spheroidal in shape for all times. The 
main result is an ODE governing the evolution of the 
droplet aspect ratio. The model is validated by exten- 
sively comparing predicted deformation with both previ- 
ous theoretical and numerical studies, and with experi- 
mental data. In particular, the experimental results by 
Ha and Yang 8] , which were obtained with extremely low 
medium conductivities are well captured by the simula- 
tion with the finite-time charge-relaxation model. The 
model is used to analyze the effects of parameters and 
stresses on the deformation characteristics. The results 
demonstrate clearly that in different regimes according to 
the sign of Qt, the stresses contribute qualitatively dif- 
ferently to deformation. Last but not least, this work lays 
the foundation for the study of a more complex problem, 
namely, vesicle electrodeformation and relaxation. This 
problem is the pursuit of our future work. 
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V. CONCLUSIONS 

In conclusion, we have developed a transient analy- 
sis to quantify droplet deformation under DC electric 
fields. The full Taylor-Melcher leaky dielectric model is 
employed where the charge relaxation time is considered 
finite. In this framework, instantaneous charge relax- 
ation is treated as a special limiting case. The droplet 



Appendix A: Appendixes 

The functions /ii(^o) ~ /i5(Co) in Eq. ([M]) are given 
in the following expressions: 



/ii(^o)= / 3^,drj, (Al) 
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/l2(eo) = ^ {/ - 3 ) ^ /l4(eo) - -C0i?3(^0) 



(A2) 



(A4) 



/i3(4o) = /iU4oj, (Acij 



/l5(eo) 



J?3(Co) G3(^o)G5 (?o) - G; (eo)G5(Co) 



/ii(eo) + eoff3(eo) 



The functions f 21(^0) — /24(Co) and F in Eq. pSal) are given in the fohowing expressions: 

, 1.2 f iv'-mv'-i) . 



drj. 



(A5) 



(A6) 



/22(eo)= €0/11(^0) 



'^3(Co) / 772 JT^ d77 + 3^oi?3Uo) 



1 



3,72 ■ 

2 — 



(A7) 



/23(eo)=eo/ll(eo) 



49(1^3gg)G3(^o)g;(go) 
307V 



^^3(^0) 



(l-3y?2)(e2_3^2^2^27y4) 



drj 



(A8) 



/24(?o)=eo'(l~eo"')^ 



3r/2 - 1 



3772 



■.dr], 



(A9) 



where 

/25(eo) = 



^ = -3(/25(^o) + /26(eo)/M, 



/22(eo) (A^."l)/l2(eo) + /l3(eo) /• 3772^1 eo 

3^0 / 7^2— ~ 72 — 



6/ll('?o) Mr/l4($o) + /l5('?o) 



,2\2 



(2g2- ^2 3^2) 



(AlO) 



■rf77, (All) 



/26(eo) 



/23(eo) (Mr-l)/i2(eo) + /i3(eo) 49(l-3^2)G;(eo) , Co 



307V 



(2g2-^2_i)(l_3^2)2 
(^2 _ ^2)2 



dr]. (A12) 
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